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I. 


Introduction 


Aerospace  structures  are  often  subjected  to  Impact  loadings.  Including 
foreign  objects  entering  jet  engines,  and  dropped  tools,  hall,  and  runway 
stones  Impinging  on  exposed  aircraft  components.  For  Impacts  by  blunt 
objects  at  moderate  velocities,  the  failure  of  the  Impacted  structure  Is 
usually  of  the  "structural"  type,  rather  than  local  perforation  or  contact 
failure.  This  report  deals  with  the  structural  type  of  transient  response 
due  to  Impact  loadings. 

Because  of  the  complicated  geometry  of  most  structures,  the  calculation 
of  transient  response  by  the  usual  finite-element  computer  codes  can  be  very 
tedious  and  time-consuming.  At  present,  there  are  no  satisfactory  methods 
for  convenient  calculation  of  the  maximum  transient  stress  or  strain;  de- 
signers need  a simple  method  for  this  purpose. 

In  this  report  we  develop  a method  to  generate  a single  design  curve 
for  a given  type  of  structure.  This  curve  gives  the  maximum  strain  In  the 
structure,  which  may  have  various  dimensions  and  material  properties,  due 
to  Impacts  of  different  values  of  mass  and  velocity.  As  an  example,  a design 
curve  for  a simply  supported  beam  under  a central  transverse  Impact  Is  pre- 
sented. It  Is  shown  that  a single  curve  of  a generalized  strain  e versus 
the  ratio  M of  beam  mass  to  Impactor  mass  can  represent  many  Impact  situa- 
tions. This  curve  may  be  generated  by  theoretical  calculations  with  different 
degrees  of  sophistication,  or  by  a limited  number  of  experiments. 
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In  the  development  of  the  design  curve,  we  have  studied  six  analytical 
models  for  treating  beam  impact  response.  Three  of  these  are  lumped-para- 
meter models  consisting  of  masses  and  springs;  the  other  three  are  dlstrlbuted- 
mass  models,  which  take  the  beam  vibration  modes  Into  consideration.  Some 
of  these  models  neglect  the  contact  force  and  treat  the  beam  with  the  attached 
impactor  mass  as  a free  transient-vibration  problem.  Other  models  Include 
the  contact  force  term,  trace  the  Impactor  motion  distinctly,  and  treat  the 
beam  as  a forced-vibration  problem.  It  will  be  shown  that  for  the  models 
neglecting  the  contact  force,  the  generalized  strain  e Is  a function  of  mass 
ratio  M only  and  Is  Independent  of  other  parameters.  For  the  two  models 
which  Include  the  contact  force,  e depends  on  M and  other  parameters,  but  the 
dependence  on  these  other  parameters  Is  weak,  and  a single  e vs.  M curve  can 
still  give  an  approximate  representation  of  all  Impact  situations. 

A series  of  Impact  experiments  has  been  conducted  using  beams  of  various 
materials  and  dimensions  and  a range  of  Impactor  masses  and  velocities.  By 
plotting  the  results  In  the  e versus  M plane,  a single  curve  can  again  be 
drawn  through  all  data  points,  with  approximately  + 30%  deviation. 

In  constructing  such  a curve  for  design  purposes.  It  Is  suggested  that 
a particular  analytical  model  be  chosen  for  calculation,  this  selection  to 
be  based  on  the  amount  and  reliability  of  the  Information  desired,  the 
particular  properties  of  the  structure  (e.g.,  damping,  plastic  behavior), 
availability  of  computing  facilities,  etc.  Alternatively,  a limited  number 
of  impact  experiments  can  be  conducted.  Regardless  of  whether  an  analytical 
calculation  or  a series  of  experiments  Is  used,  the  results  can  be  plotted 
as  an  e vs.  M curve,  which  can  then  be  applied  to  predict  Impact  behavior  for 
all  Impact  situations  within  the  considered  range  of  values  of  M. 
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The  analytical  models  and  methods  discussed  here  are  applicable  to 
all  values  of  M,  so  that  for  a given  beam,  the  Impactor  mass  may  be  large 
or  small.  The  proposed  design  curve  presented  here,  however.  Is  limited 
to  small  values  of  M,  or  large  Impactor  masses.  Some  other  nondlmenslonal 
parameters  may  have  to  be  developed  to  treat  fully.  Impacts  by  small  Impactor 
masses.  For  structures  other  than  a simply  supported  beam,  a similar  method 
can  be  used  to  construct  the  design  curve. 

II.  Analytical  Solutions  of  Beam  Impact 

There  have  been  many  analytical  derivations  for  beam  response  due  to 
transverse  central  Impact.  Most  of  these  are  presented  In  the  excellent 
book  by  Goldsmith  [1]. 

In  order  to  develop  a practical  design  curve,  we  have  performed  numerical 
calculations  by  many  of  these  existing  methods,  and  also  by  some  new  methods 
or  extensions  of  existing  methods.  In  the  following  we  shall  review  the 
basic  approach  of  each  and  present  the  numerical  results.  In  all  cases,  the 
problem  treated  Is  the  midspan  Impact  of  a simply  supported  beam. 

Table  I gives  an  Index  to  the  six  analytical  models  to  be  discussed. 

As  can  be  seen,  three  of  these  are  lumped-mass  models  and  three  are  con- 
tinuous-mass models. 

1.  One-Degree-of-Freedom,  Energy-Conserved  Model 

If  we  are  only  Interested  In  the  maximum  deflection  and  strain  In  the 
beam,  and  not  their  complete  spatial  and  temporal  distributions,  a simple 
one-degree-of-f reedom  model  with  a lumped  mass  and  spring  Is  most  convenient. 
Either  one  of  two  different  basic  assumptions  can  be  adopted;  the  first  con- 
serves the  energy  but  not  the  momentum,  the  second  just  the  opposite. 
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In  the  energy-conserved  model,  it  is  assumed  that  the  spatial  distri- 
bution of  the  beam  deflection  is  identical  to  the  static  deflection  of  the 
beam  loaded  by  a concentrated  force.  By  equating  the  original  kinetic 
energy  with  the  static  strain  energy  in  the  beam,  the  maximum  deflection  can 
be  determined.  In  other  words,  only  one  vibration  mode  is  considered,  with 
mode  shape  the  same  as  the  static  deflection  curve.  This  model  may  be  de- 
picted by  a spring  with  constant  Impacted  by  a mass  m^  at  a velocity  v, 
as  shown  in  Fig.  (1).  From  the  static  beam  deflection  formula 


6 = (PL^/48EI),  it  is  evident  that 


= 48EI/L' 


(1) 


The  maximum  deflection  of  the  spring-mass  system  with  initial  velocity  v 
is  then 


max  r *^1 


(2) 


This  equation  was  first  derived  by  Young  [2]  and  is  also  presented  in  [1]. 
We  shall  now  extend  this  approach  to  obtain  the  maximum  strain  in  the  beam. 
The  static  deflection  of  the  beam  is 

2 


w(x) 


< X < I 


(3; 


where  w^^  is  the  deflection  at  midspan.  The  maximum  strain,  which  occurs 
also  at  X = L/2,  is  then 


w 

% 

♦ V 


‘.'m 


max 


h ^ 
2 3x2 


x“L/2 


6h 

,2  ’'l 

L max 


where  we  have  replaced  the  static  w^^  by  w^ 


(4) 


Substituting  Eq.  (2)  into 


max 


(4),  and  introducing  the  definition  of  generalized  strain 

_ 2 

e * e (a  /hv) 
max 


we  obtain 


1^1 

f4M 


(5) 


(6) 


Note  that  the  generalized  strain  e is  a function  of  the  mass  ratio  M only. 
This  equation  may  be  plotted  as  an  e vs.  M curve  and  will  be  discussed 
later. 

It  is  also  Interesting  to  note  that  the  final  result  of  this  model 
is  Independent  of  the  beam  mass,  as  long  as  the  total  energy  is  conserved. 

For  instance,  considering  the  model  as  shown  in  Fig.  2,  with  the  beam  mass 
equal  to  m^  (or  any  other  value),  then,  if  the  kinetic  energy  of  the  impactor 
is  equated  with  the  initial  kinetic  energy  of  the  combined  mass  m^  + m^, 
the  combined  mass  will  have  an  initial  velocity 


""o  y (m^  + m^)  • 


(7) 


Aissuming  m^^  and  m2  are  subsequently  traveling  at  the  same  velocity  after 
initial  contact,  then  the  maximum  deflection  is  also  given  by  Eq.  (2).  It 
is  recognized  that  Eq.  (7)  violates  the  conservation  of  momentum;  if  no 
energy  is  dissipated,  the  impact  is  elastic,  and  the  impactor  will  rebound. 
It  will  be  shown  later  that  due  to  multiple  sub-impacts,  Eq.  (7)  is  still 
a good  assumption  for  calculating  final  deflection  and  strain. 
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2.  One-Degree-of-Freedoin»  Momentum-Conserved  Model 

This  model  conserves  the  initial  momentum  but  not  the  initial  kinetic 
energy.  Refer  to  Figure  2 where  the  beam  is  replaced  by  an  equivalent 

mass  (17/35)m^  and  an  equivalent  spring  constant  The  expression  for 

the  equivalent  mass  is  obtained  by  matching  the  kinetic  energy  of  the 
equivalent  mass  traveling  at  the  velocity  of  the  beam  midspan  point,  and 
the  kinetic  energy  of  the  beam,  assuming  that  the  deflection  mode  shape  is 
the  static  deflection  curve.  If  the  initial  impact  is  perfectly  inelastic, 
then  momentum  conservation  yields 

m2  V » (“2  ^ “l^^O 

or 

Vq  - v/(l  + ^ M)  (9) 

After  the  initial  impact,  the  two  masses  are  attached  and  undergo  a free 
vibration  with  initial  velocity  Vq  given  by  Eq.  (9).  The  maximum  deflection 
is  then 


”1  0 Y K, 

max  f 1 


17 

35 


M) 


(10) 


This  result  is  known  as  the  Cox  equation  [1],  [3].  We  shall  now  extend 
it  to  give  the  maximum  strain.  Following  the  approach  used  in  the  pre\  Lou^ 
case,  Eqs.  (4),  (5)  and  (10)  may  be  combined  to  yield 


E 


3 

. 4M(1  + 17M/35) 


] 


1/2 


(11) 


This  expression  gives  lower  values  of  £ as  compared  to  the  energy-conserved 
case,  Eq.  (6).  For  small  values  of  M,  Eq.  (11)  approaches  Eq.  (6). 
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3.  Two-Degrees-of-Freedom  Model 


In  the  one-degree-of-f reedom  models  considered  above,  the  contact 
force  and  relative  motion  between  the  impactor  and  the  beam  are  neglected. 
As  a result,  either  conservation  of  momentum  or  conservation  of  energy 
must  be  violated,  and  the  impactor  is  always  attached  to  the  beam,  allowing 
no  rebound  or  multiple  sub-impacts.  The  simplest  way  to  include  the  con- 
tact force  is  to  consider  the  impactor  as  a separate  rigid  body  and  account 
for  its  motion,  leading  to  a two-degrees-of-freedom  model. 

According  to  the  Hertzian  contact  principle,  the  contact  force  F is 
related  to  the  "approach",  or  relative  displacement  a,  by  the  relation 

F = k2  (12) 

which  is  nonlinear.  For  simplicity,  we  shall  approximate  this  relation  by 
a linear  one. 


The  equations  of  motion  of  this  system  (Fig*  are  then 


B - KjWj  - 0 

™2'*2  ^ ® 


> for  a > 0 


(15) 


while  the  beam  and  impactor  are  in  contact.  During  periods  when  the  two 
are  separated,  the  beam  vibrates  freely,  so  that 


17 

35 


“i"i 


^ for  a < 0 


(16) 


The  second  condition  corresponds  to  free  flight  of  the  impactor.  Using 
the  initial  conditions 


Wj^(O)  ■ W2(0)  * 0 
Wj^(O)  * 0 
W2(0)  ■ V 

we  can  obtain  a solution  of  the  above  equations,  as  detailed  in  Appendix  A. 

Having  determined  the  motion  of  the  beam  w^^(t),  the  maximum  value  of 

the  deflection  w^  may  be  substituted  into  Eqs.  (4)  and  (5),  yielding  a 
max 

corresponding  value  of  e which  is  then  plotted  as  one  point  in  the  e vs.  M 
diagram.  Fig.  5.  It  is  demonstrated  in  Appendix  A that  e,  so  obtained,  is 
a function  of  M = mj^/m2  and  K » that  its  dependence  on  the 

stiffness  ratio  K is  weak  for  realistic  values  of  and  K^.  In  other  words, 
e does  not  change  much  as  the  value  of  K is  varied.  Therefore,  only  one 
curve  is  shown  in  Fig.  5 for  the  two-degrees-of-freedom  model. 
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4.  Clebsch  Solution 


A solution  of  the  beam  Impact  problem  which  considers  the  vibrations 
of  the  distributed  mass  of  the  beam  was  first  introduced  by  Clebsch  and 
subsequently  presented  by  St.  Venant  [1],  [4].  The  impactor  mass  is  con- 
sidered as  rigid  and  attached  to  the  beam  at  the  impact  point.  An  initial 
velocity  equal  to  the  original  impactor  velocity  is  assumed  at  the  beam 
point  where  the  impactor  mass  is  attached.  The  subsequent  free  vibration 
of  the  beam  is  then  solved. 

The  Euler-Bemoulli  beam  equation  is  used 

El  + pA  — = 0 (17) 

9x  3t 

where  w is  the  deflection  of  the  beam.  The  boundary  conditions  for  impact 
at  beam  midspan  are 

2 

w(0,t)  = ^ (O.t)  ‘ ‘ 0 (18) 

dx 


*1 

T* 

/ i 

a 

*1 


3 2 

El  ^ (L/2,t)  - (l/2)m  ^ (L/2,t).  (19) 

3x  3t 

The  last  condition  equates  the  shear  force  at  midspan  of  the  beam  with 
the  inertia  force  of  the  impactor  m2»  and  implies  that  the  beam  and 
impactor  are  always  in  contact  (perfectly  inelastic  impact).  Note  that 
local  contact  deformation  is  neglected  and  that  the  impactor  is  rigid. 


Ml. 


The  Initial  conditions  are 
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Since  the  eigenvalues  are  functions  of  M only,  and  since  t*  Is  a 
dimensionless  quantity  chosen  to  maximize  a function  of  M,  it  can  be  seen 
that  e is  a function  of  M only,  or 

E ■ e (M) 

5.  McQuillen  Solution 

The  sensitivity  of  the  infinite  series  in  c,  Eq.  (23),  to  small 
changes  in  t,  makes  it  difficult  to  find  the  maximum  value  of  the  strain. 

To  facilitate  this,  McQuillen  et  al.  [5]  have  modified  the  Clebsch  solution 
by  distributing  the  impactor  mass  and  non-vanishing  initial  velocity  of 
the  beam  over  a small  distance  d on  either  side  of  the  center  of  the  beam. 
Mathematically,  the  mass  distribution  is 


“l  L 

Q(x)  = 17  + ^ {H[x  - q - d)  ] 


- H[x  - (j+  d)  ] } 

where  H is  the  unit  step  function.  The  initial  conditions  are  then 


w(x,0) 


(x,0)  - Vq{h[x  - - d) ] - H[x  - (|  + d) ]> 


where 


2dM 


+ 1 


The  solution  is  assumed  to  have  the  same  vibration  mode  shapes  as  the  original 
concentrated-mass  system,  and  the  coefficients  of  a truncated  series  are 
evaluated  using  Galerkln's  method.  The  new  solution  in  strain  is  better 
behaved  with  respect  to  time,  so  that  evaluating  the  approximate  maximum 


by  a search  method  is  more  convenient. 


]3 


I 

In  Fig.  5,  the  results  of  this  solution  are  presented  for  a fixed  | 

value  of  d/L  (=  0.02667).  Not*  that  in  this  solution  e Is  also  a function 
of  M only,  and  all  Impact  cases  are  represented  by  a single  curve.  I 


6 . Timoshenko  Solution 

As  can  be  seen  from  Table  I,  the  most  sophisticated  model  to  be  con- 
sidered Is  the  Timoshenko  solution,  which  uses  continuous  mass  distribution 
and  includes  the  contact  force  {6].  The  static  Hertz's  law  of  contact, 

Eq.  12,  is  assumed  between  the  beam  and  impactor 

F - k2 

The  displacement  of  the  midspan  point  ol  a simply  supported  beam  due 
to  a dynamic  force  F applied  at  midspan,  using  the  Euler-Bernoulll  equation. 
Is 


Wi  = 


2L 


pAir^a^  i«l,3,5  1 


rt  2 2 2 

F(t)  sin  — (t-i)dx 


(25) 


Using  Newton's  law,  the  displacement  of  the  impactor  may  be  written  as 


w = vt 

2 m„ 


ft 

0 


rt 


F dt  dt 


(26) 


The  last  three  quatlons  may  be  combined  into  a single  nonlinear 
integral  equation  in  terms  of  the  contact  force  F 
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which  may  be  solved  numerically  by  the  small- Increment  method,  as  described 


^ in  Appendix  B. 

i 

A typical  calculated  contact-force  history  is  presented  in  Fig.  6, 
showing  the  phenomenon  of  multiple  sub-impacts.  The  contact-force  pre- 
dicted by  the  two-degrees-of-freedom  mass-spring  model  for  the  same  problem 
is  included  for  comparison. 

! Due  to  the  nonlinearity  of  Eq.  (12),  the  beam  deflection  and  strain 

are  functions  of  at  least  three  parameters,  M,  k^,  and  v,  or 

e - F(M,k2,v)  . 

Thus,  a two-parameter  family  of  e vs.  M curves  must  be  constructed  to 

1 . represent  all  impact  situations.  It  can  be  shown,  however,  that  the  de- 

! pendence  of  e on  k„  and  v is  weak,  so  that  a single  curve  gives  a good 

I z 

approximation  of  all  Impact  cases.  In  order  to  demonstrate  this,  a para- 

! 

j ; metric  calculation  was  made,  with  different  values  of  the  following 

F ' 

: V parameters:  flexural  rigidity  El,  contact  stiffness  k , and  inq>act  velocity 

V.  The  actual  values  used,  as  well  as  the  calculated  values  of  e , are 

f » 

all  summarized  in  Table  II.  The  same  results  are  also  plotted  in  Fig.  5 


as  a shaded  strip  to  indicate  the  slight  spread  of  data  in  the  c vs.  M 
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An  Inspection  of  Table  II  and  Fig.  5 Indicates  that,  for  any  value 
of  M,  the  maximum  change  In  e due  to  doubling  other  parameters  Is  only 
11%.  In  particular,  the  largest  change  In  e due  to  a 100%  change  In  Impact 
velocity  alone  Is  only  4%.  Since  by  Eq.  (5),  ^ e.  it  may  be  concluded 

that  the  maximum  strain  Is  very  nearly  linearly  proportional  to  the  Impact 
velocity. 

The  effect  of  Internal  damping  In  the  beam  his  been  Included  In  the 
Timoshenko  solution  by  Hoppmann  [7].  We  have  made  calculations  with  similar 
damping  Included  and  found  no  appreciable  change  In  deflection  or  strain 
for  the  range  of  Impact  problems  considered  (low  mass  ratio  M) . 

The  effect  of  local  plastic  deformation  at  the  Impact  point  has  also 
been  considered  [8].  The  use  of  a plastic,  rather  than  Hertzian,  contact 
law,  however,  had  only  a small  effect  on  the  maximum  bending  strain.  Both 
Internal  damping  and  the  plastic  contact  law  as  applied  to  the  Timoshenko 
solution  are  discussed  In  detail  In  Appendix  B. 


III.  Experiments 


A series  of  experiments  was  conducted  to  compare  with  the  theoretical 
calculations.  of  various  dimensions  and  materials  were  Impacted  at 

midspan  by  blunt  (25.4-mm  contact  radius)  cylindrical  steel  Impactors  lireu 
from  a 25.4-mm  bore-diameter  compressed-air  gun.  The  beams  were  supported 
by  rounded  (3.2-mm  radius)  blade  supports  positioned  as  close  as  practical 
to  each  end  opposite  the  air  gun.  Four  types  of  beam  specimens  were  used: 
(1)  flat  aluminum  beams,  200  mm  span,  2.31  mm  depth,  and  three  widths, 

76.5  mm,  40  mm,  and  30  mm;  (2)  an  aluminum  beam,  152  mm  span,  20  mm  width, 
and  12.7  mm  depth;  (3)  a steel  beam,  98  mm  span,  15.1  mm  width,  and  9.0  mm 


16 


T 


depth;  (4)  a graphite-epoxy  laminated  beam,  96.5  mm  span,  38.1  mm  width,  and 
1.6  mm  depth.  The  graphite-epoxy  specimens  were  fabricated  from  0.132  mm-thlck 
Hercules  AS3501  prepreg  tape  at  a lay-up  of  (+  45/O2/+  45)g.  Bending  strain  was 
monitored  by  a 1/8- In  (3.2-mm)  strain  gage  (Micro-Measurements,  Inc.,  types 
EA-06-125AD-120  and  EA-13-125AD-120)  mounted  on  the  beam  directly  opposite  the 
Impact  point,  and  recorded  on  a Tektronix  Type  565  oscilloscope.  The  Impact 
velocity  was  determined  by  recording  two  stroboscopic  exposures  of  the  Impactor 
just  before  Impact,  measuring  the  distance  between  the  Images,  and  timing  the 
Interval  between  exposures  using  the  oscilloscope  and  a photodiode  system.  The 
relative  error  In  the  velocity  measurement  Is  estimated  at  + 1%  and  in  the  strain 
measurement  at  + 3%. 

Another  series  of  experiments  was  conducted  using  a drop-weight  apparatus 
to  Impact  the  steel  beam  described  above,  and  another  steel  beam,  197  mm  span, 
15.9  mm  width  and  depth.  Impact  velocity  was  determined  from  the  drop  height 
In  this  case. 

Figure  7 shows  two  of  the  measured  strains  as  functions  of  time.  These 
two  strain  curves  were  recorded  on  the  same  beam  Impacted  by  the  same  projectile 
at  two  different  velocities.  As  can  be  seen  from  the  figure,  the  general  shapes 
of  the  curves  are  similar,  and  the  magnitudes  are  proportional  to  the  impact 
velocity. 

The  results  of  our  experiments,  as  well  as  of  a few  experiments  conducted 
by  others  [9,10],  are  summarized  In  Table  III.  They  are  also  plotted  in  Fig.  8 
In  the  e vs.  M plane.  These  results  will  be  discussed  in  the  next  section. 

Also  conducted  were  a series  of  impact-to-failure  experiments  on  graphite- 
epoxy  beams  of  various  sizes  and  lay-ups.  These  are  discussed  in  detail  In 


Appendix  C 
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The  flexural  rigidity  El  of  each  type  of  graphite-epoxy  beam  specimen 
used  in  the  Impact  experiments  was  determined  by  a static  three-point  bending 
test.  In  each  test  the  load  and  deflection  were  measured  and  the  flexural 
rigidity  was  calculated  from  the  static  beam  flexure  formula.  The  flexural 
rigidity  was  also  calculated  from  laminated-plate  theory  by  assuming  that 
El  approximately  the  same  value  was  obtained  In  each  case. 


4 > 
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IV.  Discussion  and  Design  Application 

1.  Comparison  of  Results 

Among  the  six  analytical  solutions,  the  four  that  neglect  the  contact 
force  all  yield  relations  of  the  form  e = e(M);  that  is,  the  generalized 
strain  is  a function  of  M only.  The  other  two  solutions  show  that  e depends 
weakly  on  other  parameters.  In  Fig.  5,  the  e vs.  M curves  from  the  two 
one-degree-of-freedom  models  and  the  McQulllen  solution  are  all  distinct 
curves.  The  curve  from  the  Timoshenko  solution  is  shown  as  a narrow  strip 
indicating  the  spread  of  values  of  e due  to  changes  in  other  parameters. 

The  spread  in  e calculated  by  the  two-degrees-of-freedom  model  due  to 
variation  of  the  stiffness  ratio  K » small,  however,  that  these 

results  are  plotted  as  a simple  line.  The  Clebsch  solution  is  not  shown  in 
Fig.  5. 

The  Timoshenko  solution  has  lower  values  of  e as  compared  with  the 
McQuillen  solution.  In  the  former,  the  impactor  may  rebound,  retaining 
part  of  the  total  energy,  whereas  in  the  latter,  the  Impactor  retains  no 
kinetic  energy  when  the  maximum  deflection  and  strain  are  reached.  If  the 
Clebsch  solution  were  shown,  it  should  be  slightly  higher  than  McQuillen' s 
because  Clebsch 's  conserves  the  total  energy;  McQuillen 's  solution  balances 
the  momentum  of  the  Impactor  and  a. small  portion  of  the  beam,  thus  losing 
part  of  the  initial  kinetic  energy.  The  portion  of  the  beam  that  acquires 
initial  velocity  is  small,  thus  the  loss  of  energy  is  also  very  small. 

The  same  discussion  can  be  applied  to  the  lumped»mass  models.  The 
energy-conserved  model  gives  the  highest  strain,  and  the  momentum-con- 
served  model  gives  the  lowest  strain. 
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It  can  also  be  seen  from  Fig.  5 that  the  continuous-mass  models  give 
higher  strain  than  the  lumped-mass  models.  This  may  be  attributed  to  the 
fact  that  In  the  higher  modes  of  beam  vibration,  which  are  not  considered 
by  the  slnq>ler  models,  more  bending  strain  is  produced  per  unit  of  strain 
energy. 

The  experimentally  measured  strain  is  in  general  lower  than  the  calcu- 
lated value.  In  Fig.  9,  the  measured  and  calculated  strains  at  the  back  sur- 
face of  the  beam  at  midspan  are  compared.  The  general  shapes  of  the  curves 
are  in  agreement.  The  peak  strain  calculated  by  the  Timoshenko  solution  is 
approximately  36%  higher  than  the  measured  one.  This  discrepancy  could  be 
due  to  a combination  of  material  damping,  plastic  effects.  Inaccuracy  in 
material  properties,  and  Inaccuracy  in  strain  measurement. 

It  is  Interesting  to  note  that  within  each  solution  method,  or  following 
one  series  of  experiments,  a single  e vs.  M curve  can  be  constructed  to  rep- 
resent approximately  all  impact  situations.  This  result  implies  that  for 
values  of  M less  than  about  two,  the  most  important  parameter  in  structural 
response  is  the  mass  ratio  M;  the  contact  stress  has  direct  influence  on  the 
details  of  the  multiple  sub-impacts,  but  has  little  effect  on  the  maximum 
strain.  It  is  also  demonstrated  that  the  strain  is  directly  proportional  to 
the  impact  velocity. 

2.  Design  Applications 

Additional  research  is  needed  to  ascertain  the  reasons  for  the  dis- 
crepancy between  the  analytical  models  and  the  experiments.  None  of  the 
analytical  models  can  be  used  with  certainty  to  precisely  predict  strain. 
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The  present  results,  however,  are  already  useful  to  designers  in  predicting 
trends  and  approximate  values. 
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The  designer  may  use  the  present  curves  of  the  one-degree-of-f reedom 
models  for  beam  Impact  problems.  He  may  choose  the  Timoshenko  model,  and 
add  whatever  necessary  refinements,  such  as  plastic  contact  stress  and 
material  damping,  and  obtain  a new  c vs.  M curve.  He  may  also  choose  to 
conduct  a series  of  experiments  at  low  Impact  speeds,  each  experiment  at  a 
different  value  of  M.  The  resulting  e vs.  M curve  can  then  be  used  for  pre- 
dicting the  strain  due  to  higher-speed  impacts. 

To  use  the  design  curve  in  Figs.  5 or  8,  compute  first  the  bcam-to- 
Impactor  mass  ratio  M,  then  locate  the  corresponding  value  of  e from  the 
curve.  The  maximum  bending  strain  can  then  be  found  from  Eq.  (5).  If  failure 
of  the  structure  Is  assumed  to  be  associated  with  the  maximum  strain,  then 
an  impact-to-fallure  experiment  may  be  conducted  and  the  calculated  maximum 
strain  for  this  impact  condition  can  be  considered  as  the  ultimate  dynamic 
strain. 

The  present  paper  has  concentrated  on  transverse  impact  of  simply 
supported  beams,  but  the  same  procedure  can  be  easily  applied  to  other  types 
of  structures.  In  general,  a separate  e vs.  M curve  must  be  constructed  for 
each  given  type  of  structure.  For  instance,  the  problem  of  central  impact 
of  clamped  and  simply  supported  orthotropic  plates  is  currently  being  solved 
by  the  same  method  and  will  be  presented  in  a separate  report. 

The  analytical  methods  discussed  here  can  be  used  for  impact  problems 
with  any  value  of  M.  The  design  curve,  e vs.  M,  however,  is  a single  curve 
only  in  the  region  of  small  values  of  M,  M < 2.  For  larger  values  of  M 
(smaller  impactor  mass),  some  other  dimensionless  parameters  must  be  developed 
to  represent  the  Impact.  This  is  currently  being  pursued. 
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V. 

List  of  Symbols 

2 

a 

m 

/El/pA 

b 

m 

width  of  beam 

El 

m 

flexural  rigidity  of  beam 

F 

m 

contact  force 

h 

m 

depth  of  beam 

^^2 

= 

Hertzian  contact  stiffness 

a 

equivalent  beam  spring  constant 

K2 

m 

linearized  contact  stiffness  (spring  consta 

K - 

h'h 

■ stiffness  ratio 

L 

m 

length  of  beam 

“l 

m 

mass  of  beam 

”2 

m 

mass  of  Impactor 

m 

e 

m 

(17/35)m^  * effective  mass  of  beam 

M 

m 

t 

m 

time 

V 

s 

velocity  of  Impactor  before  Impact 

^0 

a 

initial  velocity  immediately  after  Impact 

w 

a 

w(x, t)  = beam  deflection 

(- 

w(x)  • static  beam  deflection) 

’^1 

a 

w^(t)  ~ deflection  of  beam  at  midspan 

(~  deflection  of  equivalent  beam  mass 

(<■  static  deflection  at  midspan) 

- 

maximum  midspan  deflection 

max 

”2 

a 

w^Ct)  displacement  of  Impactor 

X 

- 

beam  axial  coordinate 

z 

a 

distance  from  beam  neutral  axis 

a 

a 

contact  approach  - W2  ~ 

6 

- 

static  midspan  deflection  of  beam 
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e 

s 

e (t)  = bending  strain  at  midspan 

e 

maximum  value  of  e(t) 

max 

e 

m 

2 

generalized  strain  “ e a /hv 

max 

pA 

- 

linear  density  of  beam 

T 

- 

2 2 

dimensionless  time  “ 4ta  /L 

^1 

- 

eigenvalues 

• 

3 

time  derivative 

i 

1 

I 

f I 
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Table  III.  Beam  Impact  Experiments 


Beam  material 
and  dimension, 

L X b X h (mm) 

Impactor 

mass, 

ni2 

Mass 

ratio, 

M-mj^/m^ 

Impact 

velocity, 

v(m/sec) 

Maximum 

strain, 

E X 106 

Generalized 

strain, 

E 

.816 

0.1315 

2.44 

3100 

1.895 

aceei 

.816 

0.1315 

3.44 

4600 

1.990 

102x15.1x9.0 

1.043 

0.1029 

2.44 

3650 

2.231 

(drop  weight) 

1.043 

0.1029 

3.44 

5300 

2.293 

1.724 

0.0622 

2.44 

4200 

2.567 

Steel 

.0567 

1.893 

3.90 

1200 

0.458 

102x15.1x9.0 

.0567 

1.893 

6.40 

2000 

0.465 

.110 

0.978 

3.44 

1550 

0.670 

(air  gun) 

.110 

0.978 

5.82 

2600 

0.665 

Aluminum 

.0567 

1.93 

4.39 

1500 

0.475 

152x20x12.7 

.0567 

1.93 

6.92 

2300 

0.462 

(air  gun) 

.110 

0.998 

3.96 

1900 

0.667 

.110 

0.998 

6.19 

3050 

0.686 

Aluminum 

.110 

0.908 

3.11 

2300 

0.892 

203x76.2x2.3 
(air  gun) 

.110 

0.908 

6.34 

4380 

0.997 

Aluminum 

.0567 

0.923 

2.90 

1650 

0.822 

203x40x2.3 

.0567 

0.923 

7.01 

3900 

0.803 

(air  gun) 

.110 

0.477 

3.05 

2700 

1.278 

.110 

0.477 

6.74 

5100 

1.092 

Aluminum 

.0567 

0.692 

2.87 

1810 

0.912 

203x30x2.3 

.0567 

0.692 

7.07 

4300 

0.877 

(air  gun) 

.110 

0.357 

3.32 

3300 

1.433 

.110 

0.357 

6.80 

5980 

1.269 

Steel  [9] 

1.819 

2.16 

.408 

143 

0.506 

775x25.4x25.4 

1.819 

2.16 

1.27 

483 

0.548 

(drop  weight) 

1.819 

2. 16 

1.78 

700 

0.568 

4.314 

0.909 

.381 

278 

1.053 

4.314 

0.909 

.762 

672 

1.271 

4.314 

0.909 

1.01 

983 

1.396 

Table  III  (continued) 


Beam  material 
and  dimensions 
L X b X h (mm) 


Impactor  Mass 
mass,  m2  ratio, 

(kg) 


Impact  Maximum  Generalized 

velocity  strain,  strain, 

v(m/sec)  e x 10^  - a^ 


Steel,  Ref.  [10] 
(8-ft  section  of 

44.5 

3.05 

3.46 

680 

112-lb/yd  railway 
rail) 

108.6 

1.25 

3.46 

1300 

(drop  weight) 

213.6 

0.635 

3.46 

1720 

0.0144 

0.599 

4.75 

2400 

Graphite-epoxy 

0.0144 

0.599 

9.88 

5400 

Composite, 

0.0144 

0.599 

13.08 

6200 

AS 3501, 

0.0281 

0.307 

4.72 

4100 

9o«  3x3o*  Xxx • D 
fair  2un) 

0.0281 

0.307 

6.08 

5000 

0.0281 

0.307 

7.04 

5100 

0.0568 

0.152 

4.33 

4500 

0.0568 

0.152 

6.72 

7500 

0.110 

0.0785 

3.36 

5400 

0.110 

0.0785 

3.97 

5900 

0.110 

0.0785 

7.56 

11200 

0.795 

0.487 

2.44 

1840 

Steel, 

0.795 

0.487 

3.46 

2400 

197x15.9x15.9 

1.05 

0.368 

2.44 

2125 

(drop  weight) 

1.05 

0.368 

3.46 

3100 

1.73 

0.223 

2.44 

2575 

1.73 

0.223 

3.46 

3800 

(a)  before  impact 


(b)  at  maximum  deflection 


Figure  1.  One-degree-of- freedom  model  of  beam 
Impact,  conservation  of  energy. 


(a)  before  Impact  (b)  immediately  after  Impact 

Figure  2.  One-degree-of-f reedom  (Cox)  model  of 
beam  impact,  conservation  of  momentum. 


Approach,  a 

Figure  3.  Linear  approximation  of  Hertzian  contact  stiffness 


Figure  4.  TVo-degrees-of-freedom  model  of  beam  impact. 


Figure  5.  Generalized  strain  C,  vs.  mass  ratio  M curves  as 
calculated  by  five  different  methods  for  simply 
supported  beams  subjected  to  central  Impact. 


(a)  Timoshenko  solution 


F(KN) 


(b)  Two-degrees-of-freedom  model 


Figure  6.  Contact  force  vs.  time  as  calculated  by  two  methods. 

(El  = 224.1  N-m^,  L = 152  mm,  = 1.744x10^^  N/m^^^ 
K2  =•  7.0x10^  N/m,  m^^  = 0.1096  kg,  m^  = 0.110  kg, 

V = 6.19  m/sec). 


X 106 


V = 6.19  m/sec 


V = 3.96  m/sec 


Figure  7.  Experimentally  measured  bending  strain  at  midspan 

vs.  time  for  two  Impact  velocities.  (Aluminum  beam, 
152x20x12.7  mm;  steel  Impactor,  0.110  kg;  M • 0.998) 


Graphlte/epoxy  beam 


Figure  9.  Comparison  of  Calculated  and  experimentally  measured 
bending  strain  at  midspan  vs.  time.  (Aluminum  beam, 
152x20x12.7  mm;  steel  impactor,  0.110  kg; 

V = 6.19  m/sec;  M ■ 0.998.) 


Appendix  A - Two-Degrees-of-Freedom  Model  of  Beam  Impact 

Presented  in  this  appendix  are  the  details  of  the  solution  for  the 
motion  of  the  two-degrees-of-f reedom  model  of  beam  Impact  depicted  in  Fig.  4. 
Equations  (15)  and  (16),  which  govern  the  motion  of  this  system,  are  repeated 
here  for  convenience,  substituting  the  symbol  m^  for  the  effective  mass  of  the 
beam,  (17/35)m^. 


“e^l  ■ "2^  *^l''l 

“2^2  * ^2^'^ 2 " ■ 0 


if 

(in  contact) 


(Al) 


m w,  + K,w,  ■ 0 
el  11 


if  w > w 
ir 

(during  separation) 


We  shall  first  solve  the  former  pair  of  equations,  which  govern 
the  motion  of  the  system  when  the  spring  is  compressed,  l.e.,  when  the 
Impactor  is  in  contact  with  the  beam.  The  natural  frequencies  of  this 
system  are 

2 m^K2  + m2(Kj^  + K^)  + + m2(K^  + K^)]^  - Am^m^K^^K^ 


2m  m 
e 2 


(A2) 


The  solution  for  the  displacements  due  to  free  vibration  in  teirms  of 
the  initial  displacements  and  velocities  is 


where 


w^(t)  = ACj^  sin(o)^T  + (|!^)  + BC^  sin(u2T  + ^>2) 
W2(t)  = A sin(u^T  + (fi^^)  + B sin(u2T  + <l>2^ 


= tan 


tan 


-1  ■ W2(0)C2  ■ 

J>^(0)  - *2 (0)^2  ""j 

-1  r ''1^°^  “ W2(0)Cj^  - 

0^(0)  - W2(0)C^  “2 


Wj^(O)  - W2(0)C2 
<^l(Ci  - C2)cos 


(A3) 
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Wj^(O)  - W2(0)Cj^ 

***2 


C 


1 


Sy 


and  T * t - tQ  is  the  time  relative  to  the  instant  t ■ tQ  for  which  the 
initial  displacements  and  velocities  are  know. 

During  separation,  the  latter  pair  of  equations  in  system  (Al)  de- 
scribes the  motion  of  a simple  one-degree-of-freedom  mass-spring  system 

(m  , K, ) in  free  vibration  and  of  the  Impactor  mass  m_  in  free  flight, 
el  ^ 

The  solution  of  the  equations  is 


w (0) 

w^(t)  = — — sin  + Wj^(O)  cos  Qt 
W2(t)  = W2(0)  + T ^2(0) 

where  Q = [K, /m 
1 e 


(AA) 


The  calculation  of  Eqs.  (A3)  and  (A4)  is  accomplished  with  the  aid 
of  the  computer  according  to  the  following  procedure 
1.  The  initial  conditions 

Wj^(O)  - W2(0)  = 0 

Wj^(O)  - 0 

W2(0)  - V 

where  v is  the  Impact  velocity,  are  substituted  into  Eqs.  (A3).  These 
equations  are  solved  for  the  displacements  at  small  intervals  of  time.  The 
strain  in  the  beam  may  be  derived  from  the  displacement  W2  through  Eq.  4. 

If  and  when  w^  becomes  larger  than  W2»  separation  occurs  and  the  final  values 
of  Wj^,W2,Wj^,  and  W2  are  saved  for  the  next  step  of  the  procedure. 
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2.  The  saved  values  of  the  displacements  and  velocities  are  sub- 
stituted as  initial  conditions  into  Eqs.  (A4).  These  equations  are  also 


solved  at  small  intervals  of  time,  until  becomes  greater  than  i.e., 
a second  collision  takes  place.  The  values  of  '^2  again 

stored. 

3.  For  the  second  collision  the  stored  values  are  used  as  initial 
conditions  for  Eqs.  (A3),  which  are  solved  until  separation  again  occurs. 

4.  The  solution  then  alternates  between  Eqs.  (A3)  and  Eqs.  (A4) 
according  to  the  relative  values  of  w^  and  i.e.,  depending  on  whether 
beam  and  Impactor  are  separated  or  in  contact. 

For  large  values  of  the  mass  ratio  M (greater  than  about  1.0),  only 
one  collision  takes  place.  For  smaller  values  of  M,  several  collisions  may 
occur,  as  the  Impactor  mass  m2  follows  the  beam  mass  m^  as  it  deflects. 

A listing  and  sample  output  of  the  FORTRAN  computer  program  used  to 
calculate  the  equations  are  Included  at  the  end  of  this  appendix.  The  sample 
problem  treated  has  the  following  parameters.  ' 

m^^  = 0. 1 kg 

m2  = 1.0  kg 

- 2.0  X 10^  N/m 

K2  = 20.0  X 10^  N/m 

V = 2.0  m/sec 

L = 0.20  m 

h = 0.01  m 

Displacements,  velocities  and  contact  force  are  calculated  at  0.1  psec 
intervals  and  printed  out  every  tenth  interval. 
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It  is  Interesting  to  consider  a dimensionless  form  of  Eqs.  (Al) . 
Defining  the  quantities 


K^/K2 

M = 

mj^/m2 

T = 

o)t//k 

> 

3 

^2  ' 

0)  W2/V 

then  Eqs.  (Al)  may  be  written  as 


Observation  of  these  equations  Indicates  that,  for  small  values  of 
the  parameter  K (K<<1),  the  response  in  terms  of  maximum  deflection  or 
Wj^  of  the  beam  is  not  significantly  dependent  on  the  value  of  K.  Further, 
the  peak  bending  strain  e , which  is  assumed  proportional  to  the  maximum 
deflection  (Eq.  (4))  is  also  relatively  independent  of  K.  Therefore,  M 
is  the  only  significant  parameter  in  determining  the  response  of  the  beam 
to  impact,  for  small  values  of  K.  This  is  in  agreement  with  our  numerical 
calculation  and  experimental  results  discussed  earlier.  In  most  actual 
beam- impact  problems,  the  appropriate  value  of  K is  quite  small,  since  the 
resistance  of  the  beam  to  local  indentation  K2  is  much  larger  than  its 
deflection  stiffness 


38 


ThoOdF  ■ --  CtiMPutER  PUGRAM  for  SnlVlNG  'TwO-DEGKtE.-ur-fHrk.UuM 
MODEL  OF  bEAM  IMPACT 


BEAM  UEnOING  STIFfN 

beam  effective  mass 

ACTUAL  RFAM  MASS 
HFRT^  CONTACT  STUENESS 


2 


PEAL  M1»M^#KI#K2»K»M»L»MB  

common  r l>Y?>Vl/V9»UMEGAl#r)MEGA?»0MEGAN#ni^Y12»Y21»Y22>C/0/PMl  1 > 

1 PHI2  

PI*3. 14159265 

^UO  FDRMATClaX***TlME'’»1.0X»”Tr*»l2y»"Y2”»13X»**Vr'»t3X»”V2"»lUX»'*APPKUACM 


1**/7x»"hlL.  VEL."»7X>"k  0RCE*’>6X/”STRATN'') 
iNPUT  DATA  


Vs?.0 

MB»U.l  


K 1 =2  * UtU6 


K?»20t0E06 

1=0.2 


HeO.Ol 

INITIAL  CONDITIONS 


Y ILsO, 


NM  I 

C(-Ek,AM  = SuRT(K  1/Ml  ) 
TEPMl='l*K2+M2*(Kl+<2) 

TEHi'2  = TLPMl**2-4.n*Ml*M2*Kl*K2 
“I  rH^r?  = 5wHT(TtHM2) 

PMECjA1  = CTERM1+TLPm?)/(2.0*M1*M21 


0MEGA2  = t TEKMl-TEPvi?)/i2.o*Ml*M2  1 
nMEC.Al=SOPT(UMEGAt  ) 


OME  GA2*SQHT  I UMEUA'/ 
CmEuAh=SDKTCK2/M2) 


GnTu3Ul 
1 AIPHA-Y2-YI 
IFUlPHA.LT.O*  JGOtoS 

CALLCOinACT.(T«TZER(]#Y10*Y20#V10»V20#IFLAQ) 

“■  aiFma«T2-Y1 
PR  VEL*V2-V1 

I K (M06(i#NSKIP) *Nr.0)GUT0300 

f.oru30i 

3 c«LL  Apart  CT»tzERn»Yio» Y20/Vi0»Vf0» tFlag»mu 

Al,PnA»YJ-Yl 


JK(MLU(I#NSKIP)iNP«0)GuT 
KVo  MjKV AT(1X#E12.5#8(2X»E13 
31.1  lONTlNilE 

F X = 6.0*H*Y1/L**5 
F sK2*ALVHA 


CULUMF'  Mi-Dl'iGS 
^ hTf. t 2£)T 
?01  F f'Kl;A  I ( IHl  ) 

HHlM  i?U0 

3U2  F-l.  ilvl  HjO»’T/YI»Y2»VI/ V2» 
300  KinTTfuJe 

ppm  ao#k»h 

TTuF 

FNP 


“SMtmuTTNr  At^ART  tTTTTTRO/  Y 1 0»  Y20/  V10»V20#irtA6#MULTI  M ) 

common  Y1»Y2* Vl/tf9»Pl»P2»DMEGAN»Ytl#Yt?»Y21»Y22»C»U/PHIl»eHI2 

Sn«5»in(Bmegan*ct-tzero)  ) 

CN«C0S(0MEGAN*(T-TZER0)) 


Yl»VU)*bN/l]MLGAN  + YlO*CN 

VI* V m»CN-Y10*0MEGAN»SN 


Y?»Y20+«20*( r-TZERO/ 
V2»V20 


ALPHA-Y2-Y1 

check  WMlTnrR  btAM  AND  imRaI 

IF(  ALPHA. LT.O,  )KETIJRN 


A 

MOLT IM*MULTIM+1 


U = 
YtO«Yl 


rOR  HAVE  RE-COLLIOED 


TrRMl«Yt0-Y90*C 
TrRM9«CV10-VZ0*C J/Pi 
PHl2«ATAN2(TERMl*TrRH2) 


Y1 2»Y?2«D 

I SfSTNtPl*(T*T2ERn>»PHll) 


rn* 


Yl«Yi l*SltYl2*S2 


Y9»Y91*S1*T7Z*SZ 

C1«P1 *COS(Pf  *CT-TZFRO)fPHlf ) 

C>*P?*C0Slp>*n-T7PRa)4PHl?V 

V1*Yn*Cl+Yl2*C2 

l/9iT91*tl?rP2*C7 


0£.*U.«  •121UV6C.  ul  •Iviiu'itkL  01  •^0u61ot~04  •/2V«0sc.  uU  •i0gx^it£  U4  •37092ek~04 


f 


Appendix  B.  Timoshenko  Solution  - Numerical  Scheme,  Damping  and 
Plastic  Contact  Effects 

In  this  appendix,  the  details  of  the  solution  of  Timoshenko's  equation 
for  beam  Impact,  Eq.  (27),  are  presented.  Also  Included  are  discussions 
of  the  application  of  Internal  damping  effects  and  a plastic  contact  law 
to  the  Timoshenko  method,  as  well  as  a procedure  for  averaging  the  computed 
strain  over  a small  distance  for  more  exact  comparison  with  strain  measured 
during  Impact  experiments  using  a strain  gage. 

Numerical  Scheme 

Eq.  (27),  which  Is  repeated  here  for  convenience, 


“•  fe] 
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F dt  dt 
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2 2^  2 
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F(t)  sin  — (t-T)  dr 
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may  be  numerically  solved  by  applying  the  small- Increment  method  In  which 
the  contact  force  F Is  assumed  to  be  constant  during  any  time  Increment  At. 
Expanding  the  above  Integrals  to  compute  the  force  during  the  nth  time 
Interval  F^,  we  obtain  the  following. 
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cos  — (n-j)AT  - cos  — (n-j+l)  At 


44 


(B2) 


If  the  contact  force  is  approximated  as  a linear  continuous  function 


Computer  calculation  of  numbers  raised  to  fractional  powers  Involves 
the  use  of  the  logarithm  and  antllogarlthm  functions,  each  of  which  requires 
the  machlne-tlme-consumlng  evaluation  of  a truncated  series.  For  the  sake 
of  computer  efficiency,  Eq.  (Bl)  was  first  rewritten  as  a cubic  equation  In 
so  that  all  terms  are  raised  to  Integral  powers  only.  The  resulting 
equation  Is  solved  for  each  successive  time  Increment  by  means  of  Newton's 
Iterative  method. 

Having  found  the  contact-force  history,  one  may  determine  the  maximum 
beam  deflection  and  maximum  bending  strain,  both  of  which  occur  at  midspan. 
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Note  that  the  bending-strain  coefficients  converge  much  more  slowly 


than  the  deflection  coefficients  E 


n-J+1* 


Strain  Averaging 

In  the  beam  impact  experiments,  strain  in  the  beam  was  measured  using 
resistance-wire  strain  gages.  These  devices  measure  strain  not  at  a point, 
but  over  a finite,  though  small,  area.  In  order  to  make  a more  rigorous 
comparison  between  experimental  results  and  theoretical  calculations,  the 
computed  strain  may  be  averaged  over  a distance  corresponding  to  the  length 
of  the  strain  gage. 

The  strain  due  to  bending  for  any  point  on  the  beam  may  be  expressed  as 


c(x,t)  = [ A (t)(-l)  ^ sin 

i-1,3.5  ^ ^ 

where 

A.  (t)  * — ^ [ F(t)  sin  (t-i)dT 

m^a"^  •'0  L 

Let  c be  the  average  strain  within  a distance  6 on  either  side  of  the 
avg 

impact  point  x » L/2.  Then 


‘avg‘'>  ■ 2? 


L/2+6 
L/2- 6 


c(x,t)dx 


Substitution  and  integration  leads  to 
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(B8) 


The  average  strain  may  thus  be  calculated  from  the  contact-force  history 
determined  by  the  small- Increment  solution  using  the  following  equation 
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For  the  typical  case  of  a 1/8-lnch  (3.175  mm)  strain  gage  attached  to 
a beam  of  span  96.5  mm,  <S/L  - 0.01645.  The  following  corresponding  values 
may  be  computed  for  the  averaging  factor. 
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0.9996 

0.9889 

0.9470 

0.9028 

0.6239 

It  may  be  seen  from  the  table  and  Eq.  (B8)  that  the  effect  of 
averaging  the  strain  Is  to  diminish  the  contributions  of  the  higher  modes 
of  beam  vibration.  For  cases  where  the  Impactor  Is  massive  compared  to 
the  beam  (low  mass  ratio  M),  however,  the  beam  response  is  largely  restricted 
to  the  lowest  several  modes  of  vibration.  Thus,  averaging  the  strain  for  the 
cases  studied  made  little  difference  In  the  value  of  the  strain  at  any  given 
time  (less  than  1%  decrease) 


Internal  Damping 

Hoppmann  [7]  has  extended  the  Timoshenko  solution  of  beam  Impact  to 
Include  the  effects  of  Internal  and  external  damping  and  elastic  foundation. 
If  Internal  damping  only  Is  considered,  the  equation  for  the  contact  force 
corresponding  to  Eq.  (27)  Is 
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where  c^  Is  the  Internal  damping  coefficient. 

Application  of  the  small- Increment  Integration  method  to  this  equation 
yields  the  following  equation  In  F^,  the  force  during  the  nth  time  Interval. 
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where 


n-j+1 


I ~T-  |exp[-6  (n-j)AT] 
i=i,3,5  r [ 


1^  sin  0^(n-j)AT  + cos  0^(n-j)ATj 
- exp[-  6^(n-j+l)AT]  sin  0^(n-j+l)AT  + cos  0^(n-j+l)A-^  ^ 


where  p is  the  number  of  the  highest  underdamped  mode.  Note  that  when  damp- 
ing is  absent,  6^  = 0 and  the  above  system  reduces  to  Eqs.  (Bl)  and  (B2). 

The  last  term  in  Eq.  (Bll)  is  also  the  expression  for  the  deflection  at 


midspan. 


The  expression  for  the  bending  strain  at  midspan  is 
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Overdamped  modes  were  not  considered  for  two  reasons.  First,  the 

lowest  overdamped  mode  is  quite  high  for  a moderate  amount  of  damping. 

In  a numerical  example,  Hoppmann  uses  a value  of  the  fractional  damping 

factor  B “ 6 /p  for  the  fundamental  mode  of  B,  = 0.002;  for  this  case, 
n n n 1 

the  lowest  overdamped  mode  is  n = 23.  Even  for  Bj^  ■ 0.008,  the  lowest 
overdamped  mode  is  n * 13.  For  low  values  of  beam-to-in5>actor  mass  ratio  M, 
these  higher  modes  are  not  significantly  excited.  Second,  overdamped  modes 
are  nonosclllatory , and,  once  excited,  contribute  to  the  response  of  the 
system  for  only  a short  period  of  time. 

Figure  (Bl)  illustrates  the  effect  of  a rather  large  amount  of  internal 
damping  on  the  bending  strain  calculated  by  the  Timoshenko  solution.  Note 
that  the  general  shape  of  both  curves  (damped  and  undamped)  is  quite  similar, 
indicating  that  the  first  mode  of  beam  vibration  is  not  significantly  affected 
by  damping.  The  smaller  oscillations  due  to  the  higher  modes,  however,  are 
noticeably  attenuated  after  several  cycles. 


Plastic  Contact  Law 

The  above  formulations  of  the  Timoshenko  solution  are  based  on  Hertz's 
law  of  contact 

F = k2 

a static  relation  which  treats  the  two  bodies  as  elastic.  As  suggested  by 
Crook  [8]  and  Barnhart  and  Goldsmith  [11],  a contact  law  may  also  be  used 
which  treats  the  beam  and  Impactor  as  behaving  plastically  at  the  Impact 
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The  following  assumptions  are  made; 

1.  The  initial  elastic  behavior  before  the  onset  of  yielding  is 
neglected. 

2.  The  normal  stress,  or  flow  pressure,  is  constant  over  the  area 
of  contact  during  plastic  indentation. 

3.  The  plastic  flow  of  the  material  in  the  contact  zone  is  not 
considered. 

4.  Restitution  of  the  plastically  deformed  region  is  elastic. 


These  assumptions  lead  to  a force- indentation  law  which  may  be 
written  as 


F = k2  a 


for  a > 0,  (approaching  phase) 


(B13) 
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for  a < 0,  (rebound  phase) 
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where  ^2  = 2itR  Pq,  R is  the  radius  of  the  impactor,  p^  is  the  dynamic 

flow  pressure  of  the  beam,  F and  a are  the  values  of  the  contact  force 

mm 

and  Indentation  at  the  instant  of  maximum  penetration,  and  is  the  depth 
of  the  permanent  crater.  Defining  the  recovery  ratio  n = a^/a^,  then  the 
latter  relation  may  be  written  as 
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The  contact  law  thus  described  is  represented  in  Fig.  (B2) . In  the 
beam  impact  problems  studied,  multiple  sub-impacts  occur;  therefore, 
repeated  indentations  must  be  considered.  For  simplicity,  the  permanent 
crater  created  by  one  sub-impact  was  neglected  in  the  next  sub-impact.  How- 
ever, if  during  the  elastic  recovery  in  a particular  sub-impact  the  contact 
force  again  increases,  plastic  deformation  is  assumed  to  begin  again  (Fig.  B3) . 

This  contact  law  may  be  substituted  for  Eq.  (12)  in  the  Timoshenko 
solution  and  the  response  of  a b in  to  plastic  impact  may  be  calculated. 
Choosing  values  for  the  recovery  ratio  = 0.5  and  the  plastic  stiffness 

constant  k2  = 2.0  x 10^  Ib/in,  the  solution  of  a typical  problem  using  the 
plastic  contact  law  was  compared  with  a solution  using  the  Hertzian  elastic 
law  (Fig.  B4) . Note  that  the  shapes  of  the  two  curves  are  quite  similar  but 
that  the  effect  of  the  plastic  contact  law  is  to  reduce  the  maximum  bending 
strain. 

Computer  Programs 

Listings  and  sample  outputs  of  the  FORTRAN  computer  programs  used  to 
calculate  the  response  of  the  beam  according  to  the  above  methods  are 
presented  in  the  following  pages. 

The  first  listing  presents  the  program  SMINCl,  which  Includes  the 
calculations  for  strain  averaging  and  the  effects  of  internal  damping. 
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The  values  of  the  parameters  of  the  sample  problem  are: 

m^^  = 0.1  kg 

m2  • 0.1  kg 

E * 2.1  X 10^^  N/m^ 

L - 0.20  m ' 

b - 0.012m 

h - 0.012  m 

k2  = 2.4  X 10^^  N/m^^^ 

V “ 6.0  m/sec 

6 « 0.002  m (strain-averaging  distance) 

» 0.02  (damping  ratio) 

The  contact  force  is  calculated  at  1.0-psec  intervals  and  substituted 
into  Eqs.  (B4)  and  (B5)  (or  Eqs.  (B9)  and  (B12)  if  strain  averaging  and 
damping  are  Included)  to  obtain  the  deflection  and  strain  at  midspan. 

The  second  listing  presents  SMINC2,  which  includes  the  effect  of  plastic 
deformation  at  the  impact  point.  The  same  problem  as  described  above  is 
again  treated,  with  the  additional  parameters, 
n = 0.5 

k2  = 1.0  X 10^®  N/m  . 
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Figure  Bl.  Comparison  of  midspan  bending  strain  calculated  by 
Timoshenko  small-increment  method  with  and  without 
internal  damping.  (Aluminum  beam,  152x20x12.7  mm; 
steel  Impactor,  0.110  kg;  v = 6.19  m/sec;  M = 0.998) 


Contact  Force,  F Contact  Force 
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Figure  B3.  Complication  of  plastic  contact  law  caused  by 

second  Increase  In  contact  force  during  a single 
sub- Impact. 
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Figure  B4.  Midspan  bending  strain  calculated  by  Timoshenko 
solution  using  Hertzian  and  plastic  contact  laws. 
(Aluminum  beam,  152x20x12.7  ram;  steel  impactor, 
0.110  kg;  V “ 6.19  m/sec;  M = 0.998). 


Appendix  C - Impact- to-Fallure  Experiments  of  Composite  Beams 


In  addition  to  the  strain-measurement  impact  experiments  described  in 
the  main  text,  a series  of  tests  were  conducted  in  which  laminated  composite 
beams  were  Impacted  to  failure.  This  appendix  presents  the  details  of  these 
experiments  and  demonstrates  the  construction  of  an  impact  design  curve  for 
a particular  structure  by  using  experimental  data  alone. 


Specimens 

All  specimens  were  fabricated  from  Hercules,  Inc.,  type  AS3501  graphite- 
epoxy  prepreg  tape.  Specimen  dimensions,  lamination  lay-ups,  and  other  data 
are  summarized  in  Table  Cl.  The  flexural  rigidities  listed  in  the  table  were 
determined  from  static  three-point  bending  tests,  as  described  in  the  main  text. 

Experimental  Procedure 

In  each  Impact-to-fallure  test,  a simply  supported  graphite-epoxy  com- 
posite beam  was  repeatedly  impacted  by  the  same  mass  projected  from  an  air 
gun  at  gradually  increasing  velocities  until  failure  of  the  specimen  was 
detected  by  visual  inspection.  Impact  velocity  was  measured  using  the  photo- 
graphic technique  described  in  section  III  of  the  main  text.  Experiments  were 
performed  on  several  specimens  of  the  same  type,  l.e.,  same  dimensions  and 
lay-up. 

The  average  of  the  highest  Impact  velocity,  v.,  for  which  a particular 

a 

type  of  specimen  did  not  fall,  and  the  lowest  velocity,  v^,  for  which  any 
failure  was  detected,  is  regarded  as  the  critical  or  failure  velocity, 

^d-"^f 

V “ — 

c 


2 


(Cl) 


This  velocity  Is  used  In  characterizing  the  Impact  resistance  of  each  type 
of  specimen  for  a particular  Impact  mass  (see  Table  C2) . 


i 

t 


In  several  cases,  the  no-fallure  velocity  v^  for  a specimen  was  greater 
than  the  failure  velocity  v^  for  another  specimen  of  the  same  type;  l.e., 
some  specimens  were  more  Impact  resistant  than  other  similar  specimens. 

Such  expected  "scatter"  In  Impact  resistance  was  not  large,  however,  and  no 
special  steps  are  taken  to  account  for  It  In  the  table. 


Failure  Modes 

In  specimens  which  were  only  slightly  damaged  by  the  Impact,  the  damage 
almost  exclusively  consisted  of  rupture  of  the  beam  surface  opposite  the 
Impact  point.  This  was  accompanied  In  more  damaged  specimens  by  breakage  of 
the  Impacted  surface  and  by  delamlnatlon.  Apparently,  the  predominant  mode 
of  Initial  failure  was  tensile  fiber  breakage  on  the  surface  opposite  the 
Impact  point. 


Construction  of  Design  Curve 

The  generalized  strain  associated  with  each  of  these  failure  Impact  cases 

has  been  calculated  according  to  the  formula 

2 


'f  hv 


(C2) 


c 

where  Is  the  surface-ply  strain  at  beam  failure  for  each  specimen  type. 

We  shall  adopt  the  maximum  strain  failure  theory;  l.e.,  failure  occurs  when 
the  strain  In  the  fiber  direction  of  a ply  reaches  a critical  value,  taken 
here  as  0.0112.  Since  the  surface  piles  of  our  laminates  are  not  oriented 
at  0°,  Initial  failure  does  not  occur  on  the  surface,  but  occurs  In  the  outer- 
most 0°  piles.  The  value  0.0112  for  the  flexural  failure  strain  was  published 
by  the  material  manufacturer  (Hercules  Product  Data  Sheet  No.  832). 
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An  Impact  design  curve  may  be  constructed  by  plotting  the  experimental 
data  as  points  on  an  c vs.  M graph  (Fig.  Cl).  Note  that  these  points  are 
very  close  to  a curve  representing  the  one-degree-of -freedom  energy-conserved 
Impact  model  described  In  the  main  text.  This  may  be  a mere  coincidence. 

The  fact  that  the  experimental  points  can  be  approximately  represented 
as  a single  curve  suggests  that,  for  a particular  type  of  structure,  only 
one  Impact-to-fallure  experiment  may  be  required  to  construct  a complete  de- 
sign curve.  That  Is,  only  one  test  Is  necessary  to  fix  a point,  through  which 
the  curve  may  be  drawn  as  a straight  line  on  a log-log  plot  of  vs.  M.  The 
single  test  determines  the  value  of  the  constant  C In  the  design  curve 
equation 


e 


f 


(C3) 


Alternate  Design  Curve 

The  design  curve  of  Fig.  5 In  terms  of  e vs.  M gives  the  maximum  strain 
In  the  beam  for  a given  Impact  situation.  For  the  Impact-to-fallure  experi- 
ment, the  failure  strain  Is  assumed  known,  and  the  design  curve  of  Fig.  Cl 
tells  us  under  what  Impact  condition  this  failure  stress  will  be  reached. 

For  this  purpose.  It  Is  desirable  to  plot  Fig.  Cl  In  terms  of  a dimensionless 
Impact  velocity  v,  defined  as 


which  Is  the  reciprocal  of  e,  but  with  used  as  c.  For  a 

of  M,  the  beam  will  not  fall  at  small  values  of  v or  v.  As 

Increased  to  the  value  v , or  v , failure  will  occur,  where 

c c 

related  by 


(C4) 

given  value 

V,  or  V,  Is 

V and  V are 
c c 


V 

c 


a e . 


(C5) 
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Fig.  (C2)  shows  the  experimental  test  points  In  terms  of  v vs.  M,  which 

c 

again  can  be  approximated  by  one  curve.  This  curve  divides  the  plane  Into 
two  regions,  an  unsafe  region  above  the  curve  and  a safe  region  below.  With 
this  curve,  the  designer  can  conveniently  determine,  for  a given  Impactor 
mass  ni^,  what  maximum  velocity  v^  the  beam  can  withstand. 


Table  Cl  Impact-to-Failure  Test  Specimens 


Note:  Beam  failure  is  assumed  to  occur  in  the  outermost  0®  layer  when  its  longitudinal 
strain  reaches  0.0112. 


Figure  Cl.  Design  Curve  as  Determined  by  Experiment 


Figure  C2.  Alternate  Presentation  of  Design  Curve 
Showing  Impact-to-Failure  Experiment  Data 
for  Graphite-Epoxy  Laminated  Beams 
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Great  Lakes  Carbon  Corp. , N.Y.,  NY  10017 

(Attn;  Mr.  W.  R.  Benn,  Hgr.,  Markey  Development) 1 

Gnnasan  Aerospace  Corporation,  Bethpage,  L.I.,  NY  11714 

(Attn:  Mr.  R.  Hadeock,  Mr.  S.  Dastin).  2 

Hercules  Powder  Company,  Ine.,  CuDd>erland,  MD  21501 

(Attn:  Mr.  D.  Rug)  1 

H.  1.  Thompson  Fiber  Glass  Company,  Gardena,  CA  90249 

(Attn:  Mr.  M.  Myers) 1 

ITT  Research  Institute,  Chicago,  IL  60616 

(Attn:  Dr.  R.  Cornish) 1 

J.  P.  Stevens  & Co.,  Inc.,  N.T.,  NY  10036 

(Attn:  Mr,  H.  I.  Shulock) 1 

Kasuin  Aircraft  Corporation,  Bloomfield,  CT  06002 

(Attn:  Tech.  Library).  ......  • 1 

Lehigh  University,  Bethlehem,  PA  18015 

(Attn:  Dr.  G.  C.  Sih) 1 

Lockheed-Califomia  Co^Htny,  Burbank,  CA  91520 

(Attn:  Mr.  E.  K.  Walker,  R.  L.  Vaughn)  2 

Lockheed-Georgia  Company,  Marietta,  GA 

(Attn:  Advanced  Composites  Information  Center,  Dept.  72-14, 

Zone  42).  1 

LTV  Aerospace  Corporation,  Dallas,  TX  75222 

(Attn:  Mr.  0.  E.  Dhonau/2-53442,  C.  R.  Foreman) 2 

Martin  Company,  Baltimore,  MD  21203 

(Attn:  Mr.  J.  E.  Pamken)  1 

Materials  Sciences  Corp.,  Blue  Bell,  PA  19422  ..........  1 

McDonnell  Douglas  Corporation,  St.  Louis,  MO  63166 

(Attn:  Mr.  R.  C.  Goran,  0.  B.  McBee,  C.  Stenberg) 3 

McDonnell  Douglas  Corporation,  Long  Beach,  CA  90801 

(Attn:  H.  C.  Schjelderup,  G.  Lehman)  2 

Minnesota  Mining  and  Manufacturing  Company,  St.  Paul,  MH  55104 

(Attn:  Mr.  W.  Davis)  1 

Northrop  Aircraft  Corp.,  Noralr  Div. , Hawthorne,  CA  90250 

(Attn:  Nr.  R.  D.  Hayes , «J.  V.  Noyes,  P.  E.  I^e).  .......  3 

Roclwell  International,  Columbus,  OH  43216 

(Attn:  Mr.  0.  G.  Acker,  K.  Clayton) 2 

Rocimell  International,  Los  Angeles,  CA  90053 

(Attn:  Dr.  L.  Laekaan) 1 

Rociwell  International,  Tulsa,  OR  74151 

(Attn:  Mr.  E.  Sanders,  Nr.  J.  H.  Powell) 2 

Owens  Coming  Fiberglass,  Granville,  OH  43023 

(Attn:  Mr.  D.  Mattes).  ...........  ...  1 


75 


Mon-GovernMfit  Agencies  (Cent.) 


Rohr  Corporation,  Riverside,  CA  92503 

(Attn:  Dr.  V.  Riel  and  Mr.  R.  Elkin) 2 

I^an  Aeronautical  Conpaiqr,  San  Diego,  CA  92112 

(Attn:  Mr.  R.  Long) 1 

Sikorsky  Aircraft,  Stratford,  CT  06497 

(Attn:  Mr.  J.  Ray)  1 

Southwest  Research  Institute,  San  Antonio,  TX  78206 

(Attn:  Mr.  G.  C.  Griines) 1 

University  of  Oklahoau,  Noxnan,  OK  93069 

(Attn:  Dr.  G.  M.  Nordby)  ..............  1 

Union  Carbide  Corporation,  Cleveland,  OH  44101 

(Attn:  Dr.  H.  F.  Volk)  1 


